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We consider the finite temperature phase diagram of holographic QCD in the Veneziano limit 
(N c —> oo, Nf —y oo with Xf = Nf/N c fixed) and calculate one string-loop corrections to the 
free energy in certain approximations. Such corrections, especially due to the pion modes are 
unsuppressed in the Veneziano limit. We find that under some extra assumptions the first 
order transition following from classical gravity solutions can become second order. If stringy 
asymptotics are of a special form and there are residual interactions it may even become of 
third order. Operationally these computations imply modelling the low temperature chiral 
symmetry breaking phase with a hadron gas containing N'j massless Goldstone bosons and 
an exponential spectrum of massive hadrons. A third order transition is possible only if 
repulsive hadron interactions via the excluded volume effect are included. 
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1 Introduction 


In this paper we consider thermodynamics of QCD at finite temperature in a non-critical 
holographic model |lj and show how one can effectively compute quantum 1-loop and even 
stringy corrections by phenomenological considerations relating to the hadronic phase. We 
will then have thermodynamics at all temperatures together with an analysis of possible 
phase transitions. 

Holographic models of QCD (V-QCD) in the Veneziano limit, 


N r . 


oo 


Nf 


oo , Xf = -jJ- fixed, 


( 1 ) 


have been proposed and studied in mum, based on earlier proposals for pure Yang-Mills 
HI and the identification of the chiral condensate as the tachyon of string theory |B]. In such 
theories for a small enough number of flavours, xj < x c ~ 4 the low energy theory is QCD- 
like: it has chiral symmetry breaking, Nj — 1 massless pions (when quarks are massless) and 
confinement. 

The finite temperature study of V-QCD at zero baryon density and for massless quarks, [2] 
revealed the possibility of one or two possible phase transitions, depending on the asymptotics 
of scalar potentials: 


(a) One phase transition. In this case the theory at a specific critical temperature T c jumps 
to a black hole solution with restored chiral symmetry. The transition is a first order 
one and is at the same time a deconfinement and a chiral restoration transition. 
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(b) Two phase transitions. In this case the theory at a specific critical temperature T c 
jumps to a black hole solution which breaks chiral symmetry. The transition is first 
order and the theory is in a deconfined plasma phase with broken chiral symmetry. 
At a higher temperature there is a second transition, second order this time, and the 
system is described by a black hole solution with unbroken chiral symmetry. 

The chiral phase transition above can be characterised by a chiral condensate which is an 
exact order parameter for massless quarks. Although in the presence of quarks there is no 
order parameter for confinement, at large N c , one can use the scaling of the free energy F 
with N c . A phase in which F ~ O(N^) is a deconfined phase while a phase where F 0(1) 
is a confinement phase. In a holographic theory therefore, any phase transition in which the 
system jumps from the T = 0 saddle point (or the “thermal gas solution” as it is known) to 
a black hole saddle point, is a deconfinement transition, [j] 

In [3] a finite temperature and finite density study was made for a V-QCD model with 
two phase transitions (in case (b) above). Since then further analysis has indicated that case 
(a) above is preferred when generic properties of the mesonic radial trajectories are imposed 
mm- This is the version of the theory we will use in this work. 

As it was already stressed in |I], holographic theories in the Veneziano limit have more 
uncontrolled phenomena compared to theories in the standard ’t Hooft limit. The reason is 
that in the open string theory sector where the fundamental degrees of freedom arise (quarks) 
as endpoints of open strings, the effective coupling constant is g s Nf where g s is the closed 
string coupling g s ~ 

In the’t Hooft limit where Nf is kept fixed and of order one, g s Nf ~ jt- and contributions 
from open string loops are suppressed. This is equivalent to the fact that quark loops are 
suppressed in the QFT. In the Veneziano limit however, g s Nf ~ xj ~ 0(1) and open string 
loops are unsuppressed. This can be easily checked in a simple one loop diagram that corrects 
a propagator and pions (or in general non-singlet mesons) go around the loop: the diagram 
for a single meson is suppressed by jp but there are Nj — 1 non-singlet mesons going around 
the loop, bringing back this contribution to become of 0(1). 

In this paper we will start an investigation of such corrections. The brute-force method is to 
compute the one-loop corrections to the free-energy in the dual string theory, something that 
is in principle possible in string theory [12] . In our case however we are limited by the fact 
that we do not know the full string theory, and even if known it requires a compete solution 
at tree level in order for the full one-loop contribution to be computed. For comparison this 
is not yet known even in the best known case of holography: that of J\f = 4 super Yang-Mills. 

We will have therefore to be more modest, and we will use physical arguments to isolate 
and compute the most important contributions in the domain 0 < Xf < x c to the free energy 
of the confined saddle point: the thermal gas solution. The reason is that at tree level the 
free-energy of the thermal gas solution is 0(1) and therefore is neglected compared to the 
deconfined free energy that is O(N^). At one-loop however the correction, being proportional 
to the meson multiplicities is of O(Nj) and therefore of the same order as the deconfined free 
energy in the Veneziano limit. 

1 Of course this structure of two transitions is not specific to holography, but a similar structure emerges 
in conventional effective model approaches to QCD matter thermodynamics, see e.g. 1311 [9]. 
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The string one-loop calculation to the free energy in the confined phase can be divided into 
an infinite number of one-loop calculations where each of the string states goes around the 
loop. It depends only on the quadratic part of the tree level action and not on interactions. 
If one can diagonalize the tree-level action then the result can be given in terms of one-loop 
finite temperature diagrams of particles with given tree level masses. For the important 
meson trajectories and the present holographic theory, this was done in mam. 

For practical purposes we will split the non-singlet spectrum of the string theory in question 
(this is the one that gives the dominant contribution, the singlet spectrum contribution is 
down by a factor of NJ 2 ) as follows: 

1. The massless sector. This includes the N 2 — 1 massless pions. 

2. The massive meson sector of the four main meson trajectories. This includes the fields 
generated out of the vacuum by the three important low dimension operators in the 
meson sector: the quark mass term (massive pseudoscalars and massive scalars), the 
vector current (massive vector mesons, starting with the p) and the axial current (mas¬ 
sive axial vector mesons). The masses and widths of these mesons were computed in 
maim, where it was shown that after the lightest 2-3 mesons the rest of the masses 
are described by linear radial trajectories of the form ~ a n with a a universal 
computable constant. 

3. The full string spectrum that includes an infinite number of extra fields. Such fields 
corresponds to higher spin mesons that appear at higher levels in the string spectrum 
and therefore the lightest mass of their trajectory is higher than the four basic tra¬ 
jectories. As the detailed string theory spectrum for V-QCD is not known, we will 
parametrize such a spectrum by a density of states in order to estimate their impact 
on the thermodynamics. 


Concretely, the above is implemented as follows in this holographic setup. To begin with, 
the free energy —/ = p q (T, p = 0; N c , Nf, rn q = 0) in the QCD plasma phase is computed so 
that it is normalised to the Stefan-Boltzmann limit at T —> oo. It is valid for T m [ n < T < oo, 
were T m ; n is the temperature where the gravity solution corresponding to the QCD plasma 
phase becomes thermodynamically unstable. 

The computation of the pressure Ph{T ) down to T = 0 involves the stages 1-3 outlined 
above. First, at stage 1 Goldstone bosons with 

p,, id (T)/T 4 = (N 2 - 1) (2) 


are included. They arise as zero modes of the 1-loop computation. Comparing with p q shows 
that a 1st order deconfining transition takes place. Second, at stage 2 some states from the 
lower trajectories are included to give us ph(T). A numerical analysis shows that they only 
have a small effect and the transition still remains 1st order. Finally, at stage 3 the entire 
mass spectrum including lower Regge trajectories and with mesonic interactions is needed. 

We shall model this by an exponential Hagedorn mass spectrum jl3]-[30], see Eq. (19) 


below. Effectively, one has a mesonic string model for hadrons; these are less well developed 
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for baryons but we here address the case of zero net baryon number, p = 0. The mass 
spectrum contains a number of parameters, which are strongly constrained by how Ph{T ) 
and p q (T) connect at T c . We shall see that in this model it is easy to enforce constancy of 
both p = p/T 4 and of its logarithmic T-derivative, the interaction measure (e — 3p)/T 4 . The 
chiral transition would then be of 2nd order. However, the second logarithmic T-derivative 
thereof would then be discontinuous and the interaction measure would exhibit a sharp peak 
at T c . 

Our goal is to proceed one step further and demand also the continuity of the second log¬ 
arithmic derivative. This is actually quite complicated to achieve, as the second logarithmic 
derivatives are naturally of opposite signs. We find that to change the sign of the second 
derivative of the hadronic pressure and to accomplish required equality one has to include 
interactions in the hadron gas phase. It is enough to include repulsive Van der Waals-type 
interactions caused by the finite size of hadrons. The transition then is of third order (as 
in 2-dimensional lattice SU(1V —> oo) gauge theory [31] 32]) and actually the shape of the 
interaction measure resembles that of a Tracy-Widom distribution [55] 55] . With increasing 
Xf one sees concretely how the thermodynamics approaches that in the conformal region 
Xf > x c . 

Our analysis does not prove that the transition is of third order, just analyses a possible 
framework. Universality arguments based on the epsilon expansion (35] indicate a 1st order 
transition for Nf > 3 at fixed N c . In [5] we discuss why this conclusion might not be valid 
for the case studied here. If the transition is continuous, p{T) is analytic at all T. How the 
hadronic and plasma phases then should be connected for massive quarks is discussed in [28]. 
This is excluded in the holographic approach in which the plasma pressure p q (T ) is valid 
only for T > T m i n . Notice also that the analysis of QCD in the Veneziano limit by using 
weak coupling expansion on S' 1 x S' 3 suggests that the transition is a crossover or of very high 
order [5B] (see also m- 

One obvious extension of the above is explicit chiral symmetry breaking induced by quark 
masses. Then the massless Goldstone bosons would entirely disappear. 

Section [2] below outlines the numerical computation of p q (T) from holography [38] . and 
summarizes hadron gas formulas. Section [3] discusses a second order connection and the 
difficulties of a third order connection. The formalism of including hadronic interactions via 
the excluded volume correction is presented in Section [4] and results computed from this are 
presented in Section [5] where also the x/-dependence of the results is discussed. A simple 
modelling of what happens if quarks are massive is contained in Section [6] 

2 The pressure in high and low temperature phases 

2.1 High temperature pressure and vacuum spectrum from a holographic 
model 

To compute the pressure at temperatures above the phase transition, we apply a model for a 
5-dimensional gravity dual of QCD matter with large N c and Nf = XfN c originally proposed 
in [T] and studied in [2J and [3] (thermodynamics) and El (mass spectrum). The details of 
the model have been thoroughly exposed in M, and here we only briefly recall some of the 
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general details. 

The model builds on asymptotically AdSs metric g^, specified by ds 2 = b 2 (z)(—f(z)dt 2 + 
dx 2 + dz 2 / f(z)). In addition to gravity, the model contains bulk scalars, a dilaton A(*), a 
tachyon t(z), and the scalar potential (the zeroth component of the gauge field) 4>(z). The 
scalar 1/A sources the field theory operator TrF 2 ^. The vacuum solution A(z) is therefore 
identified with the gauge theory coupling N c g 2 (g)/(8ir 2 ). 

Furthermore, the dependence of the fields on the coordinate z in this model is constrained to 
reproduce the renormalisation group flow of the dual gauge theory in the UV (i.e., for 2 —»• 0). 
For the field A this means that 1/A(z) ~ bo log(A,z) as z —> 0, where bo = \ (11 — 2 Nf) is the 
one-loop coefficient of the beta function^} Here we also see that A is analogous to the scale 
Aqcd- F° r the tachyon r, the UV behaviour is constrained by 


t(z) 

C 


_ _LLL o 

m q z(— logAz) b o + az (— log Az) b o , 


(3) 


where 70 = I, the one-loop coefficient of the anomalous dimension of the mass operator 
in the dual field theory. Also, £ is the UV AdS radius, i.e., b(z) « Ljz as z — > 0, and 
a is proportional to the chiral condensate (qq) with a known proportionality constant j40l . 
Finally, the boundary value of the scalar potential equals the chemical potential, <£(0) = g, 
which we set to zero here. 

In the far IR, the model is required to lead to confinement. The modelling of the UV and 
IR behaviors listed above is parametrised in terms of potentials V g (X), V/o(A), k( A) and w( A) 
which appear in the action of the five-dimensional gravity coupled with the scalars discussed 
above. To determine the high temperature pressure and the vacuum spectrum, we apply the 
results of mm with the following set of potential^ 


Vg( A) = 12 


1 + “ + J 6 . 19A L + 3e- 1 /(^( 2A )4/3 v / 1 + log(1+2A) 


Vfo(A) = W 0 + A 


27 729(1 + 2A) 

24 + (11 - 2x f )W 0 


+ 


729 


24(857 - 4 6x f ) + (4619 - 1714x/ + 92x^)lU 0 


A 2 + 


120A 3 


(1 + 2A) 2 /3 


k{\) 
w{ A) 


(l+ 115 18 16a:/ A + 20A 2 ) 

VV3£ 2 

1+ 115 ~g 6a;/ A + 20A 2 ) 


2/3 


2/3 


1 + 200 log(1 + A2) • 


1 -I-log(l + A 2 ) 

200 ; 


(4) 


(5) 

( 6 ) 

(7) 


The 0(1), 0(A) and 0(A 2 ) terms are tuned so that the solutions lead to QCD beta function 
and mass anomalous dimension satisfying the standard UV expansions. The asymptotic value 

2 The constraint is imposed so that the scheme independent two-loop running is reproduced; see j2j [3j for 
details. 

3 Notice that the chosen normalization of the potentials also fixes the UV AdS radius through C 2 = 12/(12 — 

xWo). 
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Figure 1: The scaled plasma phase pressure and its first and second logT-derivatives, plotted for 
Xf = 2 and for small T (for the relations to energy density e and sound velocity c 2 , see 21 and 22). The 
potentials used are those in Eqs. 0-0. At T —>■ oo the 3 p q curve approaches 37t 2 /10 = 2.96 and the 
zero of the pressure is at T/A = 0.0925. At T — > 0 the Goldstone bosons have 3p^/(A^ 2 T 4 ) = x 2 tt 2 /30. 
Without massive hadrons there would be a 1st order transition at T = 0.125 A. Including the massive 
hadrons in 10 and 13 gives only a marginal effect in the curve marked gb. 


of Vfo in the UV, Wo, remains a parameter. Its range is 0 < Wo < 12 /xf and we have chosen 
Wo = 3/11. 

In the IR, at large A, confinement (area law) and linearity of the asymptotic glueball 
trajectories at high excitation numbers require that V g A 4 / 3 -v/log A 0. Moreover we have 
chosen the remaining parameters in 0 such that a good fit to YM thermodynamics is 
obtained, and the form of Vfo, n and w at large A (as well as the value of Wo) such that the 
phase diagram as a function of Xf is reasonable and the asymptotics of the meson spectra is 
linear with equal slopes in all sectors mm- 

To obtain thermodynamics, one searches for black hole solutions with a horizon at z = Zh- 

f(zh) = 0 , -f'(z h ) = 4vr T , s = . (8) 

Using the code in |38| thermodynamics can be computed for a given set of model functions. 
As an example, Fig. PH shows for Xf = 2 the scaled pressure p q /T A and its first and second 
derivatives with respect to logT. Note that the second derivative is negative. The pressure is 
Stefan-Boltzmann-normalised, i.e., p q /(N 2 T 4 ) — > (2 + ^ Xf)ir 2 /90 as T —> oo. It is this set of 
curves (and similar ones for Xf = 1, 2.5) that we take as the equation of state in the plasma 
phase. 

As discussed in the introduction, this holographic setting, in addition to the thermody¬ 
namics, allows also the vacuum spectrum to be determined m- The computation involves 
linearising the equations of motion around the extremal solutions of the action of the model 
and solving these linearised equations of motion. In this sense it is the lloop computation 
discussed in the Introduction. For the potentials of Eqs. 0-0 , the lowest vector masses 
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are computed to be 


2.815, 4.731, 6.076, 

0.707, 1.156, 1.477, 

0.0795, 0.1270, 0.1612, 

while the lowest scalar masses are 

2.640, 4.568, 5.857, 

0.599, 1.102, 1.405, 

0.06224, 0.1198, 0.1519 


7.188, 

(Xf = 1), 

(9) 

1.744, 

(xf = 2), 

(10) 

0.1897, 

IT 

II 

to 

(ID 


6.958, 

(x f = 1), 

(12) 

1.670, 

(x f = 2), 

(13) 

0.1802, 

(.Xf = 2.5), 

(14) 


all in units of A. The axial vector and pseudoscalar masses are less relevant since they are 
larger; the lowest axial vector masses are 4.289, 1.092, 0.1249 and the lowest pseudoscalar 
masses are 4.863, 1.173, 0.1279 for Xf = 1, 2, 2.5, respectively. Including these states and 
their radial recurrences using formulas in the following section gives the pressure marked gb 
in Fig. [l] In the relevant transition region they have a marginal effect, the temperature in 
the quark phase is so low that these massive hadronic states are hardly excited. 


2.2 The pressure at low temperature modelled as a hadron gas 

At T —> 0 the relevant degrees of freedom are the Goldstone bosons spanning the coset space 
S\J(Nf)xS\J(Nf)/SXJ(Nf). The pressure of these massless bosons is 

JV 2 T 4 30 ’ ( 

which is also shown in Fig. [lj This contribution alone, when matched with the high temper¬ 
ature contribution determined in the previous section, would lead to strong first order chiral 
transition at T = 0.125 A at Xf = 2. However, also the massive hadrons are expected to 
contribute to thermodynamics for temperatures around or higher than their masses. Above 
we saw that the calculable low-spin masses with their radial excitations have a negligible 
effect in the T range relevant for phase transitions, a complete mass spectrum, possibly also 
with mesonic interactions, is needed. 

Generally, the ideal boson gas pressure per degree of freedom is 


p h (T,p,m) 


d 3 p p 2 1 

(2t r) 3 3 E e( E ~v)/ T - 1 


T 2 m 2 

27r 2 




k =1 


m\ 


(16) 


We will set p = 0 for the rest of this paper, but we will need formulas with p to 
interactions via the excluded volume effect (see Eq. (28) below). From here one 
further 


d p h (T,m ) _ 1 Aim 3 / m \ 
8T T 4 2tt 2 ^ kT 3 1 V 1 T ) ’ 

k =1 


impose 

derives 

(17) 
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and 


d \ Ph(T,m) = J_\p 
1 dT J T 4 9^2 ^ 


oo r 4 

m 


2 t r 2 


fc=i L 


— r< , (u m \ 2 m is (u m 
T 4 °\ k T ) k T 3 1 \T 


We want to fold the pressure with the mass spectrum 

p(m , b, a , mo) = £(m) + — (e brn 9(m — mo), 

mo Vmo / 


(18) 


(19) 


where po is a dimensionless number. The degeneracy factor Nj of the m = 0 Goldstone 
bosons and the massive flavor nonsinglet states was factored out from (19). For the massive 
states with m > mo, we assume an exponential Hagedorn spectrum together with a power of 
m. Using this spectrum we get for the scaled hadron gas pressure 


Ph(T,b,a,p 0 ,m 0 ) = 


Ph 


7r 


Po 


1V 2 T 4 90 X/ + m 0 Xf J mo 


m u bm m 
dm — e 


m 


2vr 2 T 2 


/u 2 ( 


( 20 ) 


where we approximate^] the sum over k by the Hrst term k = 1 in the contribution from mas¬ 
sive states. Note that the dimensionless quantity p/T 4 can only depend on the dimensionless 
combinations T/mo and bm o- 

We will aim at matching the 1st and 2nd logarithmic derivatives of the pressure with those 
of the high temperature phase, so we compute their expressions here. By using ( |17[ ) and (18) 
we find that 


p' h ( T , b,a) = T 


d 


Ph 


dT 1V 2 T 4 


e — 3 p 
1V 2 T 4 


m 0 J , 


™ 3 


mo 


rrin 


2tt 2 T 3 


Ki 


( 21 ) 


and 


Ph( T ’ b > a ) = 


T— 

dT 


Ph 

!V 2 T 4 


(c s ~ 3)(e + p) -4(e- 3p) 


PO 2 

~ X f 
m 0 J 


m 


= - Xt / dm —- e 


bm 


m 


'mo 


m\ 


o 


2vr 2 T 3 


1V 2 T 4 
m 


m 


-K 0 [-)-2K 1 [- 


m 


( 22 ) 


Note how both of these vanish in the conformal case e = 3 p and (? s = 1/3. 

So far we have not specified the physical value of the unit of energy A. It is the same for 
thermodynamics (Fig. [l]) and the hadron masses (Eqs. (f9|) - (fT4]) ). It could be fixed in GeV 
units if one, e.g., knew T c in GeV units for Xf = 1. Its Xf dependence requires further study. 

The integrals in (]20|)-p2|), of course, blow up for T > Tnagedorn = 1 /b. However, exactly 
at T = 1/b the integrand at large m is ~ m a+3/,2+ *, where i is the order of the derivative, 
so that the ?n-integral converges for sufficiently negative a. The hadron gas pressure then 
does not diverge but approaches a constant as T —> 1/b from below. This will be the case in 
practice. 

4 Note that including only the k = 1 term provides a very good approximation: even at m = 0 the exact 
result ph(T,0) = 7t 2 /90 = 0.10966 deviates just a little from the approximate one 1 /7r 2 = 0.1013. 
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To further elaborate on the functional form of the spectral weight p(m), we note first that 
experimental hadron data cannot fix the form of the spectral weight. Simply, the available 
range of masses is too small to, e.g., separate a power from exponential. For example, over 
the range 1 < m /GeV < 2 numerically |22l [2H] 

4.52 exp(2.76 m/GeV) «---- —- 7 - exp(5.75 m/GeV ). (23) 

7 7 ((m/GeV) 2 + 0.25) 5 /2 11 v ’ 

We also note that a large number of papers have been written on p(m), starting from the 
classics by Hagedorn m Huang-Weinberg |14j and Frautschi m- Some examples are as 
follows: Ref. [16,] computes p(m) in the bag model, [17J writes “confining phase is consistent 
with effective string theory in which conformal symmetry and modular invariance play a 
significant role”, fl8] for the first time discusses a separate density for mesons and baryons, 
but believes that they should be equal to exponential accuracy, m tries to connect hadronic 
and plasma phases like here, [20]-[22] believes that the mesonic and baryonic densities should 
be different and fits both of them, [23] doubts the empirical validity of the exponential mass 
spectrum and [25] updates mass spectrum fits and includes the chemical potential in the 
hadron gas discussion. From all the work on p(m) it is obvious that there is no unique 
parametrisation for it. 


3 Connecting hadron gas with plasma 

We now keep fixed the plasma phase thermodynamics, plotted in Fig. [l] for Xf = 2. The 
solution is chirally symmetric, i.e., corresponds to the zero value of the bulk tachyon r. The 
pressure of the plasma phase vanishes at T/A = 0.0925 for Xf = 2. In the hadron phase 
the pressure is given in (20) and depends on a number of parameters. The question then is 
how QCD dynamics connects these two curves and what this implies for the properties of the 
hadron phase. 

We shall attempt to make the transition as continuous as possible. For this one needs to 
satisfy the matching conditions, in increasing order: 


• 1st order transition: only p continuous, 

• 2nd order transition: Both p and p' continuous, 

• 3rd order transition: p, p' and p" continuous. 


To begin with, in the hadronic phase at low T one at least has massless Goldstone bosons, 
which concretely arise due to chiral symmetry breaking. Their T —> 0 contribution to 3 ph, 
Xj7t 2 /30, is also plotted in Fig. PlI together with a marginal effect of the lowest calculable 
masses. The curve denoted by "gb” in the figure simply continues to the plasma line and 
gives a 1st order chiral (and deconfining) transition when one moves from one pressure curve 
to the other. The transition temperature would be at T g b = 0.125 A, somewhat higher than 
0.0925 A (where the pressure of the plasma phase goes to zero). 
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The model formally has 5 parameters: mo, b, a, po anqjthe value of T c . The minimum mass 
mo can be determined holographically by a separate computation, see Eqs. (f9|)-(fl4|) . This is 
conceivably possible also for the Hagedorn temperature 1/b, but this has not yet been done. 
Further, a is expected to be negative and T c has to be somewhat above the point at which 
p q vanishes in Fig. [lj Even though the number of parameters is in principle high enough to 
match all conditions, it is not clear that a solution exists. In fact, we shall see that a third 
order solution is found only if interactions in the hadron gas are taken into account. This 
happens a la Van der Waals by including the effects of finite size of hadrons. 

Even if we ultimately determine the parameters a, po and T c by numerical fitting, it is 
useful to have a simple toy model to estimate their values. To obtain the functional form of 
Eq. (19), consider the eigenvalues N of the operator X^=i n -^W- The degeneracy of 

the eigenvalue is (d — 2 is the number of transverse dimensions) ( 12 ] 


P(N,d) = 


1 fd- 


V2V 24 


[d- 1)/4 


N' 


d+1 

4 exp 


2vr 


d- 2 


6 


-N 


(24) 


We change this to mass density by P(N, d)dN = p(m)dm, N = a'm 2 : 
p{m) = V2Va i (^^j 


d— 1 
4 


a'm) 2 exp 2\ 


Id-2 


nv a'm 


(25) 


Choosing d = 5 and mo = 1/y a' this is of the form (19) with 


v/2 

po = —- ~ 0.177, a = — 2 , bniQ = 7 tv 2 ~ 4.44 . 
8 

Our final 3rd order numbers will not agree with this toy model. 


(26) 


3.1 Second order transition 

A second order transition, i.e., a transition where p and p' are continuous, is very easy to 


obtain. As an example, consider the model (26). Then one has two quantities to determine: 


the minimum mass mo and the transition temperature T c . One first determines mo = mo(T) 
numerically by requiring continuity of p'. Then inserting this to the continuity condition 
for p. one determines the valuc^of T c = 0.5987, giving finally m q = mo(T c ) = 3.158. The 
thermodynamics so obtained is plotted in Fig. [2] One observes that the minimum mass 
obtained from the thermal fit is close to the directly determined one in ([9]) but somewhat 
bigger. 

To have an idea of the range of acceptable parameter values, take the computed value 
mo = 2.8 from ([9]). Assume a has some fixed value. Then we have 3 parameters to de¬ 
termine, T c , b, po, but only two equations, equality of p/T 4 and its logT-derivative. By 


5 With the understanding that also the continuity of pressure is a constraint for the parameters. Then also 
the critical temperature, i.e., the temperature where the constraints are evaluated, is a free parameter. 

6 The numerical values of dimensionful quantities here and below are given in units of A, unless stated 
otherwise. 
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Figure 2: Left: Thermodynamics at Xf = 1 with the toy model parameters in (26) and the fitted 
values T c = 0.5987, mo = m m ; n = 3.158. Right: The 2nd derivatives as functions of T. Note the 
opposite signs, corresponding to the sharp peak in the interaction measure in the left panel. 


T c (p 0 ;a) b (Po;a) 




Figure 3: Fitted values of T c (p 0 ;a) (left) and 6(po;a) (right) for Xf = 1 and mo = 2.8. 


demanding that the derivatives coincide one first determines b = 6(T c ,po). The equation 
Ph(T c , b(T c , p 0 ),Po) = Pq(T c ) then gives T c = T c (p 0 ) and finally b = b(T c (p 0 ),p 0 ) = b(p 0 ). The 
outcomes for T c and b are plotted in Fig.[3]as functions of the normalisation of mass spectrum 
po and for various fixed values of a. 

3.2 Third order transition with pointlike mesons 

The next stage is to also require that the 2nd derivatives be equal when we move from ph 
to p q . Since pressure in the plasma phase is taken to be fixed and p" q < 0, one also has to 
change the sign of p ^ from positive to negative, see the right hand plot in Fig [ 2 } If the second 
derivatives are equal and negative at some T c , p'! must have a zero somewhere. At the same 
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Figure 4: The second derivative in (22) with mo = po = Xf = 1, evaluated for a = —10 and a = — 6 
and (curves from bottom) b = 0.5, 1, 1.5, 2. For b > 2 the negative part disappears. The curves end 
at T = 1/6. The zero of pfc(T, 6, a), if any, is at T « 0.5 largely independent of the values of 6, a. 


point the interaction measure p' h will have an extremum and the sharp peak has disappeared. 
To see the sign of p'l(T, b, a) in (22) one can write it in the form 


pl(T, b,a) = 


Po 


jia+l 


2vr 2 rrin +1 


Imo/T 


dyy a+ 3 e bT V [yKo(y) - 2/T 1 (y)] . 


(27) 


The sign is determined by the combination yKo(y) — 2K\(y) which has a zero at y = 2.3864. 
At small y this combination is negative and behaves approximately as — 2/y — 2y log(y), while 
at large y it is positive and behaves approximately as e~ y \Jity/2. Thus the lower limit ttiq/T 
allows only positive values if T/m o < 1/2.3864 = 0.419. To have a desired negative value the 
integral must probe the negative region below y = 2.3864 by having T > 0.419 mo- Not much 
of it is needed as seen from the numerical plots in Fig. [4] (where mo = 1), p/(T, b, a) < 0 
for T > 0.5 almost independent of the values of b and a. At this temperature then first 
derivative, i.e., scaled e — 3 p has a maximum and the sharp peak in the interaction measure 
has disappeared. 

For Xf = 1 one has mo ~ 2.8 and one can expect p^ to become negative only for T > 
0.5 mo ~ 1.4. This would push the maximum of the interaction measure to values of T much 
larger than those encountered earlier, which are about 0.6. It is physically obvious that the 
hadron gas phase cannot be thermodynamically stable up to such large T. One concludes 
that with the present hadron gas model, pointlike hadrons, the sign of p'l cannot be changed, 
the cusp in the interaction measure cannot be removed and the transition cannot be made 
of higher than 2nd order. 


4 Third order transition with excluded volume corrections 

Let us then check if including hadron interactions via the excluded volume correction |26l 
E3 [28] would make it possible to make the transition of third order and to get rid of the 
cusp in the interaction measure. One replaces V —> V — vqN where vq = 1/Tq is the volume 
of a single meson and N the number of mesons. What is the new pressure p(T,p), given 
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the pointlike boson pressure po(T,p ) in (16)? Here p is the chemical potential associated 
with the meson number IV; i.e. n = N/V = dp/dp. The chemical potential associated with 
conserved baryon number is taken to be zero. 


According to a related model in [25] the effective hadron volume is / 7rr 3 = (7.93/GeV) 3 = 
1 /Tq so that To = 0.126 GeV. In QCD T c ~ 0.15 GeV while in our units T c j A ~ 0.5 (where 
we reinstated the unit of energy A). Thus A ~ 2 T c ~ 0.3 GeV so that an expected magnitude 
at Xf = 1 is To/A ~ 0.42. We shall find below by matching the pressures of the two phases 
for Xf = 1 that T c /Tq ~ 3, To/A « 0.25. 

As shown in detail in 123, Section II, p(T, p) is obtainecQas a solution of the transcendental 
equation 

p{T, P) = Po{T, p - P(T, p)) • (28) 


Taking the partial derivative with respect to p of 


n 0 (T,p- p(T, p)) 

n(T,p) = 


1 + n 0 (T, p - ^j p(T, p)) 


one obtains for the number density 

(29) 


Thus, to apply this, one has to solve p(T,p) from (28). 

In the Maxwell-Boltzmann (MB) approximation eP( E ~0 ;$> 1, and only the k = 1 term in 
the series (16) contributes, 


p 0 (T,p) = )^ /T = Tn 0 (T,p ) 

The p dependence is a simple exponential and p8]) becomes 


P(T, p) 
Po(T, p) 


= exp 


Po(T,p ) p{T,p) 
TTq po(T,p ) 


= exp 


f n o(T, p) p(T,p) 
T 3 po(T,p) 




(30) 


(31) 


In our case p = 0, and inserting this in (31) we see that all quantities are just functions of 
T. The equation is of the general form q = e~ aq which is trivial to solve numerically. More 
formally, the solution is 

q{a) = -W (a) (32) 

a 

where W(a) is the Lambert’s function, ProductLog in Mathematica parlance. At small a 

q(a) = -W(a ) = f](-l) fc ^ + ^ ' a k = 1 - a + | a 2 - § a 3 + • • • , (33) 


k =0 


and at large a 


q(a ) = -VL(a) = - (logo — log log a + 0 (1)). 
a a 


(34) 


' In this section, po is the pressure of pointlike hadron gas, p is the pressure when excluded volume effects 
are included. 


13 














Figure 5: Left: The second derivative ( Tdr) 2 p/T 4 of the excluded volume solution (35) for to = 3 


and two values of the meson volume parameter To. The pointlike curve (18) (with k = 1) is also 
shown. Right: The p(m)-integral of the second derivative of the excluded volume solution ( |35| of (31) 
with po = Xf = 1, a = —10, Too = rn m i n = 3, T 0 = 0.5, and (curves from bottom) b = 0.5, 2/3, 1. 


Thus the pressure after excluded volume corrections is 

p(T)^p(T,m,T 0 ) = TT^w(^^j =TT^w(^ 

At small T/m, no(T) ~ exp(— m/T) is small and 


mrT /m 
2it 2 T3 K2 \T 


(35) 


p(T)= P0 (T)(l- r ^p- 


T 3 

J o 


and at large T 


p(T) = TTq 


T_ 

^0 


+-J 

(36) 

I + H- 

(37) 


It is also illuminating to evaluate the effect of volume exclusion on the number density 


at large T. In (29) there is an additional suppression factor exp(— p/(TT $)) = 0(1)Tq/T 3 , 
where we used (31) and ( f37[ ) to obtain the estimate. Therefore the excluded volume density 
at large T is simply n = 0(1 )Tq, the hadrons are densely packed but do not overlap. Related 
to this, in the excluded volume model the effective chemical potential —p(T, 0)/Tq is always 
negative and one does not get to the Bose-Einstein condensation domain n —> m where meson 
wave functions overlap. 

As a first step towards understanding the effects of excluded volume, the left panel of Fig. [5] 
shows numerical results for the 2nd logT derivative of the scaled pressure solution (135]) of 


(28) for m = 3. There the red curve corresponds to the pointlike pressure, the k = l term 
in (18). As discussed above, the 2nd derivative is negative for T > m/ 2.386 = 1.257. The 


two excluded volume curves correspond to meson volume parameter values To = 0.2 and 0.5. 
Here 0.5 is chosen so that the hadron volume is of the order of l/T 3 with T c as in previous 
figures. The small T behavior is always like in the pointlike case, the large T behavior 


shows the positive 2nd derivative following from (37) and in between there is a negative 
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region. Increasing To would extend the negative region to the right so that it asymptotically 
approaches the pointlike curve. In any case meson finite size enhances the negative region. 

As a second step, the right panel of Fig. [5] shows numerical results for the 2nd logT 
derivative of the scaled pressure solution (|35l) of (28) at Tq = 0.5 but now integrated over the 


exponential mass spectrum (19) p(m,b,a,mo). A large negative value of a is needed, here 
a = —10. Since we have an exponential mass spectrum, the computation can only be valid 
up to T = 1/6. Since plasma extends down to T = 0.5, we certainly must have 6 < 2. 

Comparing with Fig. [4] one sees that indeed excluded volume based interactions increase 
the magnitude of the negative 2nd derivative even when integrated over exponential mass 
spectrum. The parameter po can be used to increase the magnitude further. 


Above in (30) we used the Maxwell-Boltzmann approximation which, due to its explicit 
p, dependence, led to simplified computations. Without invoking the MB approximation, 
instead of q = e ~ aq , one has to solve the equation 


q = 


2_a 


^oo 1 


kaq K 2 (k S 


ETw k 2 (*sr 


(38) 


for q = q(a,m/T). One can numerically check that this improvement has only a marginal 
effect. The physics reason for this is that, as discussed above, in the excluded volume model 
the meson wave packets are densely packed but do not overlap, and one is never close to 
Bose-Einstein condensation. 

Another calculable interaction term is that among Goldstone bosons. Including terms of 
order 2 and 4 in the chiral Lagrangian and to 3 loops [39l 123] one has 


P S b(T) = W^T 4 ( 1 + 


NjT' 




144/, 


4 log^ + 0(T b 


(39) 


where the scale A p depends on the higher order couplings of the bosons. For the interaction 
measure we find 

e gb ~ 3p gb _ vr 2 2 NjT 4 ^ Ap 1\ 

N'jT* 30 Xf 108// V ° S T 4 ) ' 



In the Gell-Mann-Oakes-Renner relation (for a theory with all quark masses equal) m 2 f 2 = 
~WJ m q (qq) one has (qq) ~ NfN c so that f n scales as // ~ N c . Thus in the interaction 
measure above Nj /// ~ and it seems that this term will be dominated by effects from 
the massive states. As a side remark, it is also repulsive for T < A p , it increases p g b there. 


5 Fit of parameters and their Xf dependence 


We shall now analyse the QCD equation of state (EoS) at Xf = 1, 2, 2.5 by requiring conti¬ 
nuity of the logarithmic derivatives (T<9/<9T) n (p/T 4 ), with n = 0,1,2. We use p q computed 
from holography and the hadronic phase EoS 


p h (T,b,a,p 0 ,m 0 ) = 



+ 


„ roc 

— x / [ 
mo 1 Jm a 


dm 


m a 

mg 


e bm 


T[ 

T 3 


W 


m 2 T 
2t r 2 T 0 3 


I<2 
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Figure 6: Pressure and interaction measure with a 3rd order phase transition at Xf = 1, i.e., with 
Nj = N% massless Goldstone bosons. Stable phases are continuous, metastable ones dashed. The 
maximum of interaction measure is at T = 0.672 in the hadron gas phase, hadron gas is the stable 
phase for T < 0.771 and ends at T = 1.0. 



Figure 7: Second logT derivative of p/T 4 and the sound velocity squared with a 3rd order phase 
transition at Xf = 1. One sees concretely how the 2nd derivative is continuous at T c = 0.771 but the 
third derivative jumps. The sound velocity squared is continuous at T c but not its derivative. 


where W(a) is as discussed above in (32). The first and second logT derivatives can be 
computed analytically but lead to lengthy expressions. 

Of particular interest is to see what happens at larger Xf, for values approaching the 
lower limit of the conformal region. The naive argument comparing the number of degrees 
of freedom in the conformal limits at T = 0 and T = oo gives the estimate x c = 4. The 
precise value for the potentials of the bulk action used here can be obtained numerically as 
explained in [Tf and is x c = 3.187. The value 2.5 is already rather close to this, but not 
yet in the Miransky scaling region x c — Xf <C 1, where the bound state masses and critical 
temperatures are expected to decrease as ~ exp (— const /^/x c — Xf). 

The outcome of a numerical application of the excluded volume model for Xf = 1 is shown 
in Figs. [6] and [7J One plots p, its first and second derivatives and the sound velocity, related 
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Figure 8: As Fig. [6] but for Xf = 2. Hadron gas ends at T = 0.22, is the stable phase for T < 0.173, 
the maximum of interaction measure is at T = 0.147. Note that if we fix the critical temperature to the 
QCD value T c w 0.15 GeV, we have here A « 1 GeV, i.e., the numerical values for the temperatures 
are close to their physical values measured in units of GeV. 



Figure 9: As Fig. [7] but for x/ = 2. The 2nd derivative is continuous at T c = 0.173 but the third 
derivative jumps. The sound velocity squared is continuous at T c but its derivative has a discontinuity. 


to the second derivative as in Eq. (22). Thus indeed a connection between hadron gas in the 
plasma phase with a third order transition can be established. Not surprisingly, for a wide 
range of T, 0.7 < T < 0.9, the pressures of the hadron and plasma phases are very close 


to each other. The role of the repulsive interactions in the hadron gas phase was to bend 
down the 2nd derivative to negative values. It vanishes when the interaction measure has a 
maximum. The sound velocity approaches the conformal limit l/\/3 both when T —> 0 and 
T —> oo. With broken chiral symmetry and massive Goldstone bosons, c s —> 0 at T —> 0. 

The outcome at xj = 2 is shown in Figs. [8] and [9] and that for Xf = 2.5 in Fig. |10| (where 
only p and the sound velocity are plotted, 1st and 2nd derivatives are qualitatively as in 
Figs. [8] and [9]). The fit parameters with mass dimension are summarised in Table [lj Here 
mo = 777 m i n is the smallest vector mass, T(p q = 0) is the temperature at which the plasma 
pressure vanishes (this is the transition temperature for pure Yang-Mills theory), T(p q = p g o) 
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Figure 10: Left: Fit with a 3rd order phase transition for Xf = 2.5: the 2nd derivative of p is 
continuous at T c = 0.025 but the third derivative jumps. The maximum of the interaction measure 
(not shown) is a T max = 0.0209. Right: The sound velocity squared. It is continuous at T c but not 
its derivative. 

is the temperature at which the plasma pressure is p g b = -/V|7 t 2 T 4 /90, at T max the interaction 
measure has a maximum, at T c chiral symmetry is restored (the plasma phase becomes the 
stable phase), Tq is the meson volume parameter (vq = 1/Tq) and 1/6 is the Hagedorn 
temperature (the metastable hadron gas phase ends there). 

Although mo varies by more than a factor 35, the expectation that parameters with mass 
dimension scale with mo is born out to a 10% accuracy. An exception is the meson size 
parameter Tq which decreases by more than expected by mass scaling. Comparing with T c 
we have T c /Tq = 3.1, 4.3, 4.6 for xj = 1, 2, 2.5, respectively. One may note that even at 
Xf = 1 the value of To = 0.25 was smaller that the expected value 0.42. Thus with increasing 
Xf the mesons have to appear effectively larger, the interactions stronger, to bend the hadron 
gas EoS nearly continuously to the plasma one. 

The numerical value of the mass exponent in p(m ) ~ e bm can be understood by noting 
that physically the hadron gas as a metastable phase cannot be expected to extend far 
into the plasma phase. The end point is the Hagedorn temperature Th = 1/6 so that 
1/6 > T c ~ 0.3mo, in agreement with the observation 6mo ~ 3. 

The dimensionless parameters vary on a similar level: (a, po) = (—9.61, 6.25) for Xf = 1 is 
changed to (—8.5, 5.67) for Xf = 2 and to (—8.1, 6.27) for Xf = 2.5. It is not excluded that 
constant Xf independent values could be found for these. Anyway it seems that the exponent 
of the powerlike mass dependence is stably close to —8 and the weight of the massive part of 
the mass spectrum is about 6. 

Note the small range of variation in c 2 , which gets monotonically smaller when Xf increases: 
one is approaching the conformal region where everywhere c s = l/\/3. 

One may convert the above temperatures to GeV units by, for example, demanding that 
T c (xf = 1) = 0.15 GeV and by assuming that A is independent of Xf. Then A = 0.194 GeV 
and T c = 150, 34, 4.9 MeV at Xf = 1, 2, 2.5, respectively. For the determination of scales 
at large Nf and estimates of the N / dependence of the critical temperature by using other 
methods, see, for example, muMsa. 
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x f 

m 0 

T{p q = 0) 

T{p q =p g b) 

T max 

T c 

To 

b 

1 

2.82 

0.480 

0.542 

0.672 

0.771 

0.25 

1 

2 

0.707 

0.0925 

0.125 

0.147 

0.173 

0.04 

4.5 

2.5 

0.0795 

0.0124 

0.0179 

0.0209 

0.025 

0.0054 

30.3 


x f 

m 0 

o 

II 

cy 

II 

Oq 

cr 

T max 

f c 

T 0 

bmo 

1 

2.82 

0.170 

0.194 

0.238 

0.273 

0.0887 

2.8 

2 

0.707 

0.131 

0.177 

0.208 

0.244 

0.0566 

3.2 

2.5 

0.0795 

0.156 

0.225 

0.263 

0.314 

0.0679 

3.0 


Table 1: Top: Dependence of mass scales on Xf. Here mo is the minimum vector mass. All 
quantities are in units of A. Bottom: Dependence of mass scales on if scaled with m®: 
T = T/m 0 . 


6 Second order transition with massive Goldstone bosons 


If quark masses are non-zero, they break explicitly chiral symmetry and Goldstone bosons 
are massive. The hadronic pressure will then vanish at T = 0. To see what this implies 


quantitatively, take the hadron mass spectrum in (19), use the toy model parameters in (26) 
and replace 5(m) —>■ 6(m — m g b) with m g b = 0.5 ~ T c . The argument here is that for physical 
QCD T c is close to m n . Also the holographic plasma part will, in principle, change; there 
are only tachyonic chiral symmetry breaking solutions even though the quark mass is very 
small. We do not yet have these solutions available and stick to the m q = 0 plasma curves. 
The thermodynamics computed with these assumptions is shown in Fig. m 

At small T the effect, of course, is striking. One observes further that the minimum mass 
obtained from thermodynamics is now almost equal to the one from direct computation in 
(19) and that T c is only 10% above the temperature at which p = 0. Due to this closeness 


the derivatives are also larger and the peak in interaction measure even sharper. 


7 Conclusions 

We have in this paper shown concretely how a high T, p = 0, QCD plasma equation of state, 
computed from holography at vanishing quark mass, can be connected with a low T hadron 
gas phase with Nj Goldstone bosons and massive mesons obeying a Hagedorn spectrum with 
a minimum mass. In the language of holography, the leading holographic computation has 
been improved by quantum lloop and stringy corrections. The holographic computation 
implies that there is a minimum temperature for the plasma phase and, accordingly, a phase 
transition is needed. It is very simple to connect the phases with a first or second order 
transition, and we have shown how a more continuous third order transition can be achieved. 
The motivation for this is that for physical non-zero quark masses lattice Monte Carlo results 
suggest that the transition is continuous. 

Quantitatively a determining role in the hadron gas phase is played by the minimum mass 
of the hadrons (mesons in our case). We emphasize that here these minimum masses have 
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Figure 

bosons, 


11: Left: Thermodynamics with toy model parameters in (261 but with massive 
m gb = 0.5 ~ T c . The fitted values are T c = 0.5348, mo = 2.754. 


Goldstone 


also been computed from the same holographic model at T = 0 and at various values of 
Xf = Nf/N c . The entire spectrum obviously cannot be computed. We have seen here how 
it and its interactions are constrained by the requirement of as continuous a transition as 
possible. 

An outcome of the computations here is that the thermal parameters with mass dimension 
scale with the minimum mass when Xf is varied. This is expected to also be true in the 
Miransky scaling region, i.e., very close to the start of the conformal region at Xf = x c ~ 4. 

Our computation, of course, does not prove that the transition is of third order, it just 
indicates what phenomena are encountered if this is the case. Conceivably one could also 
impose continuity of the third derivative. Completely analytic expression |28lj would, however, 
require that the plasma EoS extend to T = 0 and the hadron gas one to T = oo, which is 
not possible in our model. 

Ultimately, of course, these theoretical ideas should be tested by numerical lattice Monte 
Carlo simulations, say, at N c = 3 and Nf = 3, 6, 9,.. as approximations to N c , Nf —> oo. 
Much work has been devoted to this at T = 0. The figures above should give a good idea of 
what can be expected to happen to thermodynamics when xj is increased. The overriding 
difficulty is the imposition of vanishing or small quark mass. This holographic computation 
can also be extended to non-zero /r. 
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